Chirality-induced Dynamic Kohn Anomalies in Graphene 
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We develop a theory for the renormalization of the phonon energy dispersion in graphene due 
to the combined effects of both Coulomb and electron-phonon (e-ph) interactions. We obtain the 
renormalized phonon energy spectrum by an exact analytic derivation of the phonon self-energy, 
finding three distinct Kohn anomalies (KAs) at the phonon wavevector q = u)/v, 2&f ± lo/v for LO 
phonons and one at q — cu/v for TO phonons. The presence of these new KAs in graphene, in 
contrast to the usual KA q — 2/cf in ordinary metals, originates from the dynamical screening of 
e-ph interaction (with a concomitant breakdown of the Born-Oppenheimer approximation) and the 
peculiar chirality of the graphene e-ph coupling. 
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Graphene, a two-dimensional (2D) sheet of graphite 
comprising a planar hexagonal lattice of carbon atoms, 
has attracted a great deal of recent interest with enor- 
mous amount of experimental and theoretical activity 
since its discovery. Of fundamental significance is the 
study of the interaction between electrons and phonons 
in this novel system. In particular, the many-body renor- 
malization of the phonon spectrum due to the electron- 
phonon (e-ph) interaction is an important open question. 
Experimentally, Raman scattering offers a means to mea- 
sure the long-wavelength phonon energy whereas X-ray 
and neutron scattering provide an avenue to probe the 
entire phonon energy dispersion. Due to free carrier- 
induced many-body interactions, one expects the ob- 
served phonon energy in a doped system to be differ- 
ent from the "bare" phonon energy in an undoped sys- 
tem. Physically, electrons respond to the dynamical lat- 
tice vibrations by screening the lattice potential, chang- 
ing the elastic constants of these vibrational modes and 
thereby renormalizing the phonon energy. Recently, a 
number of Raman scattering experiments [1] on extrin- 
sic graphene (where the electron density can be tuned 
by gating) have emerged, with the observed density de- 
pendence of the Raman shift (for the long-wavelength 
G-band optical phonons at the T point) pointing to the 
inapplicability of the Born-Oppenheimer approximation 
(BA) [2| in graphene. This density dependence has been 
addressed for the long-wavelength phonons (i.e. q = 0) 
using perturbation theory [J y| and density-functional 
theory |4|. 

An interesting question remains as to whether the 
Kohn anomaly (KA) [H|, which appears as a cusp in the 
phonon energy dispersion at q = 2kp for ordinary 2D 
metals, will be substantially modified in graphene due 
to its quasi-relativistic chiral band structure. The occur- 
rence of K A is entirely a many-body effect as it originates 
from the screening of the e-ph interaction by electrons, 
and correpsonds to the singularities of the phonon self- 
energy or its derivatives as a function of q. For graphite, 



the KA was recently studied in Ref. Q as a function of 
q at zero doping. Ref. [4j addresses the phonon energy 
renormalization in doped graphene, but only at q = 
relevant to Raman scattering experiments. Tradition- 
ally, KA is probed experimentally with X-ray or neu- 
tron scattering spectroscopy by measuring the phonon 
energy as a function of q, for which a systematic theory 
for the phonon energy renormalization in graphene for 
finite wavevector q > is still lacking. 

In this Letter, we present a theory for the renormaliza- 
tion of the phonon energy dispersion in graphene due to 
both Coulomb and e-ph interaction effects. From our the- 
ory, we obtain two major new results: (1) We find that 
direct Coulomb and phonon-mediated electron-electron 
interactions decouple to all orders of perturbation theory 
within the random-phase approximation (RPA), and the 
electronic collective plasmon mode does not contribute 
to the phonon energy renormalization; (2) we obtain the 
renormalized phonon energy dispersion as a function of 
q, predicting the occurrence of three distinct KAs at the 
phonon wavevector q — lu/v, 2kp ±uj/v for the LO mode 
and one at q = lu/v for the TO mode, which arise from 
the chiral structure of the graphene e-ph coupling. The 
novel pecularity that these KAs do not occur at q = 2fcp 
(as in usual metals) originates from the fact that the 
phonon dynamics cannot be neglected in the screened e- 
ph interaction, indicating the inapplicability of the BA. 

Graphene, behaving as a 2D zero-gap semiconductor, 
has a G-band optical phonon energy uiq — 200meV at 
the r point which is comparable to the Fermi energy 
Ef ~ 110 — 370meV at the usual extrinsic carrier den- 
sity n = 10 12 — 10 13 cm -2 . This approximate equality of 
the phonon and electron energy scales loq ~ ep implies 
a breakdown of the static approximation for the phonon 
degree of freedom 0] , naturally explaining the violation 
of the BA in the recent Raman scattering experiments 
[l[ since the BA implicitly assumes the phonon dynam- 
ics to be much slower than the electron dynamics. In 
addition, plasmon-phonon coupling, which has been ex- 



tensively studied in doped semiconductor systems (e.g., 
GaAs, SiC), is expected to occur whenever the phonon 
dynamics is at a comparable time scale as the electron 
motion. It follows that one cannot take the screening 
of the e-ph interaction to be simply static while keep- 
ing the screening of the Coulomb interaction to be dy- 
namic; instead, one has to take into account the dynami- 
cal screening of both Coulomb and e-ph interactions, i.e., 
direct electron-electron (e-e) and e-ph interactions must 
be treated on an equal footing Q. 

Near the Brillouin zone corner K (i.e. the Dirac point), 
graphene is described by the effective chiral Hamiltonian 
H = vcr k, where v w 10 6 ms _1 is the quasiparticle veloc- 
ity and er is the set of Pauli matrices describing the two 
A and B sublattice degrees of freedom. The quasiparti- 
cle energy dispersion as obtained from this Hamiltonian 
is £fcA = Aefc, where = vk and A is the chirality label 
representing the conduction band (A = 1) and valence 
band (A = —1); the corresponding eigenstate is denoted 
as |feA). The e-ph interaction vertex for optical phonons 
is given by gM, where g is the coupling constant 

characterizing the magnitude of the e-ph interaction, and 
M is the off-diagonal matrix 



M{q) = 





Mba^ 







(1) 



with Mab = —1 or i and Mba = 1 or i for LO or TO 
phonons, respectively, and <p q = tan -1 (q y / q x ) the azith- 
muthal angle of the momentum q. The chiral structure 
of the e-ph coupling is peculiar to graphene, describing 
the bond stretching and bending between neighbouring 
carbon atoms of the A-sublattice and B-sublattice. 

In the presence of direct e-e interaction and e- 
ph interaction, the phonon Green function is renor- 
malized by both, as represented in the diagram- 
matic language in Fig. [TJ Summation of these dia- 
grams leads to the renormalized phonon Green func- 
tion D(q,u) = D°(q,u)/[1- D°(q,u)WP(q,u)], where 
D° (q,w) = 2wo / (w 2 — Wq + i0 + ) is the bare phonon Green 
function with uq the optical phonon energy, II pp (q,ii;) is 
the RPA-screened phonon self-energy given by the fol- 
lowing equation [Fig. QJb)]: 

ipp(g, w ) = nP p ( q ,c)+nP c ( g ,c)K ! c c (g^)n^(q,c), (2) 

where V c %{q,uj) = V q /[1 - V q U c c (a, to)} is the usual RPA- 
screened Coulomb interaction [ll| (here IIq c is the elec- 
tronic polarizability) , II PP (g, u>) the bare phonon self- 
energy, and n pc (q,oj), HQ P (q,uj) are the bare "hybrid" 
bubbles with one Coulomb interaction vertex and one 
e-ph interaction vertex. In regular metals or doped 
polar semiconductors, the phonon self-energy and the 
hydrid bubbles are the same (up to a factor given by 
the e-ph coupling constant) as the electronic polarizabil- 
ity; the Dyson equation for the renormalized phonon 
Green function (Fig. [1} simply reduces to the usual 
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FIG. 1: (Color online) (a) Dyson equation for the renormal- 
ization of the phonon Green function. The zigzag lines denote 
the phonon Green function and the crosses denote the e-ph in- 
teraction vertex. The shaded bubble with two cross vertices 
stands for the renormalized phonon self-energy, (b) Equa- 
tion for the renormalized phonon self-energy. The unshaded 
bubble with two cross vertices denotes the bare phonon self- 
energy, the two bubbles with one cross vertex and one dot 
vertex are "hybrid bubbles" with one e-ph interaction vertex 
and one Coulomb interaction vertex. The wavy line stands 
for the usual RPA-screened Coulomb interaction. 



RPA series. In graphene, due to the presence of a chi- 
ral structure of the e-ph interaction Eq. JT]), the bub- 



bles rig , IIq , n o p , and Hq c are not equal to one an- 
other. In particular, IIq C and IIq P incorporate the cou- 
pling effects of the Coulomb and e-ph interactions within 
the RPA, and describe the renormalization effect of the 
phonon energy by the Coulomb interaction. Interest- 
ingly, we find that with the graphene e-ph interaction 
Eq. (JTJ) , these hybrid bubbles vanish identically: IT pc = 
Hq P = kvTgY, ikn , k Gl x (ik n )Gl +qX ,{ik n + i ?n )<k + 
qA'|M(q)|kA)(kA|k + qA') = 0, where G a kX (ik n ) = 
l/(ik n — is the quasiparticle Green function with 
S,k\ = f'k\ — £f the quasiparticle energy rendered from 
the Fermi level. Eq. then implies that Il pp is simply 
given by the bare phonon self-energy IT PP , and Coulomb 
interaction does not contribute to the screening of the e- 
ph interaction within the RPA. The phonon energy dis- 
persion, which is given by the pole of the real part of the 
renormalized phonon Green function, 



2woRen pp (g,w), 



(3) 



is therefore only renormalized by the e-ph interaction 
but not by Coulomb interaction. It follows that cou- 
pled plasmon-phonon modes do not arise in graphene, 
with the phonon and plasmon modes having separate 
branches of energy dispersion in the lu — q phase space 
despite comparable energy scales for the phonon and 
electron dynamics. In addition, direct Coulomb and 
phonon-mediated e-e interactions are simply additive, 
with the RPA-screened total e-e interaction given by 
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V C T = V q /[1 - V q ILf] + V c f/[l - D°IP Q % where K ph 
is the unscreened phonon- mediated e-e interaction [101 ] . 
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FIG. 2: Dynamic phonon self-energy n^ p (a;,u) versus x — 
q/kF at u — u>o/£f for the LO (solid line) and TO (dashed) 
modes at a density n = 10 13 cm -2 . Inset: Static phonon self- 
energy n pp (:E,w) versus x — q/kF at U = 0. 

The phonon self-energy Il pp = II PP can then be derived 
from the e-ph interaction vertex Eq. [T]as: 



ip p (<^n) = 4^ 2 EE nF(a : ) r nF f fc+gv) 



AA' k 



[1 t- AA'cos 



*<2Vi + £fcA — Ck+qX' 
'k+q + <f>h)] , 



(4) 



where the sign — (+) corresponds to LO(TO) phonons, 
the factor of 4 counts the spin and valley de- 
generacies, and cf>k = t&n~ 1 (ky/k x ) is the azith- 
muthal angle of the momentum k measured from 
q. The real part of the phonon self-energy com- 
prises two contributions ReII pp = EE P + n pp (drop- 
ping the notation 'Re' for the real part there- 
after), where n pp = 2g 2 Y J xT. k n F{S,k lJ ,)Ux{k,q) with 
fixx(k, q) — [It Acos(0 fe+g + <j> k )\ /[u + nik+ - £,k+ q x] - 
[1 T Acos(0 fe + 4>k- q )] I [w - A*£fc+ + ik- q x\- The contri- 
bution EE P is due to the extrinsic conduction band elec- 
trons whereas II PP is solely due to the intrinsic valence 
band electrons. In the renormalization of the phonon 
energy in extrinsic graphene, II pp should be subtracted 
from the total II PP to avoid overcounting of the intrinsic 
contribution, since the "bare" phonon energy (i.e. when 
graphene is undoped) , by definition, already includes the 



effect of n pp [!,l9|. Therefore, only II PP should be taken 
into account in the phonon energy renormalization by 
free carriers in extrinsic graphene. We also note that, 
for II PP , the largest energy scale is the Fermi energy ef 
whereas for Il pp it is much higher, of the order of the 
cutoff energy A = 2nhv/a for the graphene linear band 
dispersion (a = 2.46A is the graphene lattice spacing). 
Therefore, although the BA is inapplicable for extrinsic 
graphene with ujq/ef ~ for intrinsic graphene the 

BA is strictly valid since ojq/A <C 1. 



For clarity, we express our results in the follow- 
ing dimensionless quantities: wavevector x = q/kp, 
energy u = uj/ef, and phonon self-energy n pp = 
II pp /(2(7c G £F/7r) with g% c — g 2 A/H 2 v 2 the dimensionless 
phonon-mediated e-e coupling constant (^4 is the sample 
area). We have obtained the following asymptotic re- 
sults for the phonon self-energy in the limit q/kp <C 1: 
nj&,TO+(s.«) = (l/2)[l + (u/4)ln|(2-t0/(2 + u)|] + 
AII P q TO+ (x,w), where the next order correction to the 
long wavelength result is given by, 



An pp + (x, u ) 
An pp 0+ (x, u ) 



8 - 4m 2 - u 4 1 
2u 2 (u 2 -4) 2 ~ fa n 



2 — u 



24 - 12u z 



2u 2 (u 2 - 4) 2 



1 

8u' 



In 



2 + u 
2 -it 



x 2 (5) 



(6) 

At x = 0, the long wavelength result [l], 0, EH is the 
same for LO and TO phonons, the two phonon modes 
being degenerate at the T point. At finite wavevector 
x, this degeneracy is lifted with the leading-order cor- 
rection going as x 2 given by Eqs. (J5])-©. We have also 
obtained the following analytic results for the phonon 
self-energy in the static case u = 0: n P [ P :)+ (a;, 0) = 0, and 

x) T (l/4)[(2/x)Vx 2 -4 - 



n£ p o+(^o) 



-(tt/8)x0(2 



xia.iT 1 {2/y/x 2 ~4)]6{x~2). 

The above static phonon self-energy results are de- 
picted in the inset of Fig. [21 which clearly shows the 
presence of a non-analyticity at q — 2kp corresponding 
to the KA for LO phonons. We also note that this non- 
analyticity is entirely absent in the static electronic polar- 
izability [l3( of graphene. This is in contrast to the situa- 
tion for regular materials with a parabolic energy disper- 
sion where the phonon self-energy and the polarizability 
are equal up to a factor given by the e-ph coupling. This 
distinctive difference between the phonon self-energy and 
the polarizability in graphene is a direct result of the pres- 
ence of a chiral structure in the graphene e-ph coupling 
Eq. (JXJ) , which leads to a different Berry phase depen- 
dence in the expression of the phonon self-energy Eq. ([JJ 
compared with the polarizability [l3( • Therefore, KAs in 
graphene originate entirely from the special chiral struc- 
ture of the e-ph coupling Eq. ((T|). 

We have evaluated the full expression for n pp with gen- 



eral q and ui dependence analytically [14|, which is how- 
ever too cumbersome to be shown here. The main plot of 
Fig. [2] shows the evaluated Il pp as a function of x at the 
phonon energy u> = loo, from which three cusps occurring 
at vq = u>q, vq = 2ef ± too for the LO mode are clearly 
discernible. For the TO mode, there is a divergence of 
n pp at vq = ujq. These non-analyticities correspond to 
the values of q where the denominator of the integrand of 
n pp vanishes at the Fermi surface k = kp, i.e. the zeros 
of the equation lu ± vkp ± v\kp ± q\ = 0. 

With the calculated phonon self-energy, the renormal- 
ized phonon energy spectrum can be obtained by self- 
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consistently solving Eq. Q for lo, which is shown in 
Fig. [3] for the LO mode and Fig. g] for the TO mode. 
For LO phonons, three KAs which correspond to the 
non-analyticities of n pp are evident, occurring at the 
wavevector vq = loq, vq — 2ef ± Wq. For TO phonons, 
the divergence of n pp at vq = loq is removed due to the 
self-consistency condition for lo in Eq. ([3]), but the KA 
remains as a sharp but finite peak at vq = loq. In addi- 
tion, we note the KA at vq — loq for both the LO and 
TO modes, unlike the other two KAs for the LO mode, 
is independent of electron density, an interesting conse- 
quence of the quasi-relativistic linear dispersion peculiar 
to graphene. For both LO and TO phonons, our re- 
sults suggest that the phonon energy first increases (i.e., 
phonon hardening) with density up to a certain phonon 
wavevector, and then decreases (i.e., phonon softening) 
with density. The critical wavevector for this transition 
from phonon hardening to softening is different for LO 
and TO phonons, and we find numerically q ~ 5 x 10 8 m _1 
for LO phonons and q = loq/v = 3 x 10 8 m _1 (i.e., the 
KA) for TO phonons. 
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FIG. 3: Renormalized LO phonon energy spectrum co p h versus 
q at different electron densities n. The bare phonon energy 
is shown as the horizontal solid line. The range of phonon 
wavevector q shown corresponds to [0, 0.08(27r/a)] away from 
the F point. 



In the following, we provide a schematic picture for 
understanding the occurrence of the KAs. The inset of 
Fig.|¥]shows the lo — q phase space in which single-particle 
excitation can occur through virtual phonon exchange. 
The behavior of the phonon self-energy is characterized 
by six different regions of the lo — q phase space where n pp 
is analytically continuous, separated by the boundaries 
(indicated by the solid lines) corresponding to the set 
of values of (lo, q) where LT PP is non-analytic. Necessary 
(but by no means sufficient) conditions for the KAs to 
occur are given by the intersection points between the 
phonon dispersion line lo = loq and the boundaries for the 
different regions. Within the BA, phonons are treated as 
static with lo = in the expression of the phonon self- 
energy, the only intersection points therefore occur at q = 
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FIG. 4: Renormalized TO phonon energy spectrum aj p h ver- 
sus q. The legends are the same as in Fig. [3] Inset: Different 
regions for the analytical behavior of n^ p . The solid lines in- 
dicate the boundaries for these regions, and the dashed line 
shows the phonon energy lo — loq at a density n = 10 13 cm -2 . 



and q — 2kp, as in the case of, e.g., a usual metal. In 
graphene where the BA is invalid, phonon dynamics must 
be included with lo = loq in the phonon self-energy, the 
intersection points now occur at vq = loq, vq — 2ep — loq, 
and vq — 2ep + lo . For the LO mode, we find that 
KAs occur at all three intersection points. This however 
does not apply for the TO mode, and we find two equal 
but opposite contributions in the expression for n P : P - )+ 
which cancel the effects of two KAs, yielding only one 
KA at vq — loq in this case. In ordinary metals, the KA 
q = 2k f corresponds to backscattering with maximum 
phonon wavevector q Q; this KA becomes q = 2kp^fLOo/v 
in graphene because of the phonon non-adiabaticity, with 
q = 2kp — loq/v, 2kp + loq/v corresponding, respectively, 
to emission and absorption of a phonon through electron 
backscattering. The KA q — loq/v does not have an 
analogue in ordinary metals, and corresponds to forward 
scattering of the electron through absorption of a phonon. 

In conclusion, we have developed a theory for the 
interaction-induced phonon renormalization in graphene, 
and discovered new and multiple KAs in the renormalized 
phonon dispersion. The peculiarity and the distinction of 
these KAs from the usual KAs in metals are signatures 
of the renormalized dynamically screened e-ph interac- 
tion and the special chiral structure of the e-ph coupling 
in graphene. The graphene phonon energy dispersion 
can be measured with double resonance Raman scatter- 
ing or electron energy loss spectroscopy on a monolayer 
graphene, and the experimental verification of our pre- 
dictions would establish that graphene has a very unique 
e-ph many-body coupling. 

This work is supported by US-ONR, NSF-NRI, and 
SWAN SRC. 
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